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INTRODUCTION 

Field of the Invention 

The field of the invention is computational design and analysis of arrays of genomic 
data and products. 

Background of the Invention 

With the completion of the Human Genome Project, genetic research is now being 
directed towards understanding complex multigenic diseases e.g. cancer, cardiac diseases. 
Microarray technology has proven to be really useful in studying the expression pattern of 
thousands of genes simultaneously. Also with the availability of the entire genome, many 
5" tools have now been developed to generate inferences and predictions based on the entire 
genome, such as POMPOUS (Fondon et al, PNAS, 95(13)7514-9, 1998) which looks for 
potentially polymorphic genes. Efforts like Program For Genome Application (PGA) are now 
being undertaken to study hundreds of genes associated with particular diseases or 
phenotypes. As a result, researchers frequently need to compile large lists of genes associated 

10 with certain diseases, phenotypes, keywords and their synonyms. The selection of array 
elements for large gene collections typically involves: finding possible gene candidates, 
generally done using a series of keyword searches on different databases; assembling these 
several different lists obtained from various databases and trying to eliminate redundancies; 
and annotating all the genes (on the tentative list) in detail so that the researcher can know as 

1 5 much as possible about the gene. 
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The NCBI website provides a keyword search engine for various databases like 
GenBank, UniGene and LocusLink; however, the keyword search has to be done separately 
on each database. The list then needs to be combined and more importantly, the sequence 
redundancy needs to be eliminated. Eliminating the redundancy manually is not an easy task 
5 since each database has its own unique identifier. This is primarily done based on the 

researcher's experience and not all sequence redundancies are eliminated, especially for a 
large collection of genes. Additionally, the annotation for all the candidates on the list is not 
available in one place, so the researcher has to look up individual genes - a very laborious 
and time-consuming task. 

10 Websites like Genecards (Rebhan, M et al, Bioinformatics 14(8)656-64, 1998) 

(http://nciarray.nci.nih.gov/cards/) provide a database of human genes, their products and 
their involvement in diseases. However, Genecards only offers information about the 
functions of all human genes that have an approved symbol, and a few selected others. Again 
this information can only be accessed one gene at a time, and the annotation cannot be 

1 5 downloaded in any useful format for working with a large gene collection. DRAGON 

(Bouton CM et al, Bioinformatics 16(11)1038-9, 2000)(http://207.123.190.10/dragon.htm) 
lets the researcher do a keyword search on multiple databases at one time, but the output is a 
list of accession numbers and definitions in text format, which is not linked to any of its 
annotations. The tool does not let the researcher select entries from the keyword search. It 

20 does not allow moving between pages and merge lists obtained from different keyword 
searches. As a result DRAGON does not help in systematically compiling a large gene 
collection. Further, DRAGON does not include important databases like GenBank and 
LocusLink that are the most commonly used databases for searching candidate genes. None 
of these tools helps in eliminating sequence redundancies within the lists. Databases like 

25 LocusLink and Genecards attempt to integrate the unique characteristics from various 

databases and provide a broad summary on a single gene basis. Nevertheless they do not help 
in annotating a large gene collection. There is a need for a tool that comprehensively gathers 
annotation related to all these elements in one place. The annotation tool of DRAGON only 
combines information from UniGene, Swissprot, Pfam and KEGG pathway database with 17 

30 fields of annotation. However these fields do not include important fields like repeat, SNP, 
pathways, clones, etc. which would be of great value. Additionally including a number 



2 



UTSD:668 



(expression data for microarrays, purity of repeats for polymorphism) in the final annotation 
table would make it convenient for the user to extract information from the table. With more 
and more gene collections, it is also required to combine several collections of genes, 
obtained from different sources. 
5 The production of DNA microarrays can be divided into four stages: a. Selection of 

array elements and design of the probe DNA; b. Preparation of the probe DNA; c. Preparation 
of a suitable design substrate to spot the probes on; d. Deposition of array elements. The 
selection of array elements for microarrays involves assembling a large gene collection. It 
would be very valuable if the same tool (to compile a large gene collection) could be used to 

10 further design primers, look for commercially available clones (expression microarrays) and 
design resequencing probes (resequencing microarrays). Once the genes are spotted on the 
microarray and hybridized to fluorescent labeled probes, there are a number of software 
programs that help in conversion of the fluorescence of the scanned image to numbers, using 
complex mathematical corrections to extract signal from background noise, e.g. 

15 Genepix(http://www.axon.com/GN_GenePixSoftware.html) and ArrayVision 

(http://imaging.brocku.ca/products/Arrayvision.htm). These numbers indicate level of 
expression. Other programs such as GeneSpring (Silva et al, HMS Beagle: The BioMedNet 
Magazine Issue 82, 2000), Cluster Treeview (Eisen MB et al, Proc Natl Acad Sci U S A 95) 
and Spotfire (http://www.spotfire.com), help in the analysis by clustering the data together 

20 using various methods based on K-means, hierarchal or self-organizing maps. Clustering 

algorithms use the expression level data to group the various elements on the array. It would 
also be very useful to view the elements of the array with their complete annotation and 
overlay the expression level data on top of it. The data could further be selectively viewed by 
sorting on various annotation fields and the expression level data. This approach could be 

25 useful to view any large gene collection in general. With the increasing number of microarray 
experiments, it would be valuable to compare elements between different microarrays 
considering that fragments of the same gene might be represented by different sequence 
identifiers. For example, two different accession numbers might belong to the same UniGene 
cluster, representing the same gene. An artifact sometimes observed in the results obtained 

30 from an expression profiling microarray experiment is that some sequences might hybridize 
to other sequences to which they are significantly similar. This leads to false positive results 
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after a microarray experiment. Although Human Cot DNA is often used to prevent 
non-specific hybridization by blocking simple repetitive elements in genomic DNA, as shown 
in experiments to study cross-hybridization, Human Cot DNA is not very effective in 
preventing cross hybridization. ARROGANT computationally estimates the amount of cross 
5 hybridization for each sequence and tags potential genes as possible candidates for cross 
hybridization. 

Several computational tools and databases are available which may be used in the 
development of the code for working with large gene collections. Some of them are discussed 
here in brief. 

10 1. PRIMO: PRTJVIO (Li et al, Genomics 40(3) 476-85,1997) is a code that was 

developed to design primers for large-scale DNA sequencing projects. PRIMO designs 
primers (short sequences typically 20 bases long), which are used to amplify sequences (0.4 
KB- 2 KB) using PCR. PRTMO can be made to design primers to amplify a specific region. 
PRIMO can be run in batch mode and the region for the design of primers for each sequence 

1 5 can be specified separately. The parameters file (including parameters like oligo length, 
melting temperatures etc.) can be altered. The code is written in ANSI C and is available 
locally on a HP/UX computer. The code has been successfully used to design primers for the 
past couple of years and is available on the web at http://atlas.swmed.edu. This makes 
PRTMO a very important tool to design primers to amplify a large number of sequences 

20 simultaneously. 

2. BLAST: BLAST(Basic Local Alignment Search Tool) is an alignment tool to 
search for similar sequences (protein or DNA) developed by NCBI (Altschul et al, Journal of 
Molecular Biology 215(3)4-3-10,1990). It is available at 

http://www.ncbi.nlm.nih.gov/BLAST/. ARROGANT uses the BLAST output to estimate 
25 cross-hybridization for microarrays. Each element on the array is BLASTed against the entire 
UniGene database and the BLAST output is parsed to detect 65 contiguous hydrogen bond 
overlaps, used as a threshold for cross-hybridization. 

3. Rep-X: Rep-X (Wren et al, American Journal of Human Genetics 67(2)345-56, 
2000) uses the UniGene database and generates a list of repeats, hairpin and palindrome 

30 sequences. This code runs on HP/UX computer. The output of Rep-X is incorporated into 
ARROGANT to look for repeats, hairpins and palindrome sequences. 
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4. NCBI Databases: NCBI provides databases used by ARROGANT (downloaded and 
implemented locally) to annotate gene collections and find potential candidates associated 
with keywords. The databases include: a. GenBank (Benson DA et al, Nucleic Acids Res 
28(1)15-18, 2000): An annotated collection of all publicly available DNA sequences provided 
by NIH; b. UniGene (Schuler, J Mol Med 75(10)694-8, 1997): Partitions GenBank EST 
sequences into a non-redundant set of gene oriented clusters; c. LocusLink (Pruitt et al, 
Nucleic Acids Res 29(1)137-40, 2001): Integrates and provides a single query interface to 
cluster sequences and makes available descriptive information about genetic loci; d. 
HomoloGene (Zhang et al, J. Comp. Biol. 7(1-2)203-14, 2000): The database of calculated 
orthologs and homologs between all UniGene clusters by each pair of organisms. 

5. KEGG Databases: KEGG (Kyoto Encyclopedia of Genes and Genomes) (Kanehisa, 
M., Oxford University Press 2000) provides genome and pathway databases for a large 
number of organisms. ARROGANT uses (downloaded and implemented locally) these 
databases to look for potential gene candidates, their pathways and to annotate gene 
collections. 

6. Clone Databases: Commercially available clone databases include the IMAGE (G. 
Lennon et al, Genomics 33(1)151-2, 1996) Consortium, which shares high quality arrayed 
cDNA libraries and provides sequence, map, and expression data on the clones in these arrays 
to the public domain; vendors include Research Genetics, Incyte Genomics, etc. 

SUMMARY OF THE INVENTION 

The invention relates to computer-based systems and corresponding methods for the 
design and analysis of biopolymer sequence arrays. 

In a first principal embodiment, the invention provides a computer-based system for 
creating a targeted collection of sequences from a dataset comprising sequence identifiers 
corresponding to natural complex biopolymer sequences and linked to corresponding 
annotations, the system comprising: 

a) a search function which searches the annotations of the dataset according to a user- 
defined criterion and outputs a first subset of the dataset restricted by the criterion; 

b) a redundancy reducing function which compares the first subset with a first 
database correlating the sequence identifiers of the first subset with syngeneic biopolymers 
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and outputs a second subset of the dataset having reduced unique, natural complex 
biopolymer redundancy relative to the first subset; 

c) a selection function which applies to the second subset a user-defined selection 
parameter and outputs a third subset restricted relative to the second subset by the parameter; 

5 and 

d) a tabulation function which creates and outputs the targeted collection of sequences 
in the form of a data table comprising, configurable by and sortable by the sequence 
identifiers of the third subset. 

The system may optionally incorporate one or more of the following limitations: 
1 0 the criterion is selected from the group consisting of a keyword and a concept; 

the criterion is one of a plurality of user-defined criteria, and the search function 
searches the annotations of the dataset according to the criteria and outputs a first subset of 
the dataset restricted by the criteria; 

the criterion is one of a plurality of user-defined criteria, and the search function 
1 5 searches the annotations of the dataset according to the criteria and outputs a first subset of 
the dataset restricted by the criteria, wherein the criteria include multiple keywords; 

the dataset is selected from the group consisting of GenBank, Medline and KEGG; 
the dataset is one of a plurality of datasets, and the search function searches the 
annotations of the datasets according to the user-defined criterion and outputs a first subset of 
20 the datasets restricted by the criterion; 

the database is selected from the group consisting of UniGene and LocusLink; 
the database is one of a plurality of databases correlating the sequence identifiers of 
the first subset with syngeneic biopolymers, and the redundancy reducing function compares 
the first subset with the databases and outputs the second subset of the dataset; 
25 the parameter is selected from the group consisting of source, species, author and 

pathway; 

the parameter is one of a plurality of user-defined selection parameters, and the 
selection function applies to the second subset the parameters and outputs the third subset 
restricted relative to the second subset by the parameters; 
30 the redundancy reducing function outputs a second subset of the dataset which 

eliminates unique, natural complex biopolymer redundancy relative to the first subset; and 
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the system further comprises an expansion function which searches a second database 
for synonyms of the sequence identifiers of the first, second or third subset. 

In a second principal embodiment, the invention provides a computer-based system 
for creating a targeted collection of sequences from a plurality of datasets comprising 
5 sequence identifiers corresponding to natural complex biopolymer sequences, the system 
comprising: 

a) a merge and redundancy reducing function which compares the datasets with a 
database correlating the sequence identifiers with syngeneic biopolymers and creates a subset 
of the sum of the datasets having reduced unique, natural complex biopolymer redundancy 

10 relative to the sum; and 

b) a tabulation function which creates and outputs the targeted collection of sequences 
in the form of a data table comprising, configurable by and sortable by the sequence 
identifiers of the subset. 

The system may optionally incorporate one or more of the following limitations: 
15 the merge and redundancy reducing function further comprises a selection function 

which applies a user-defined selection parameter whereby the subset is restricted relative to 
the sum of the datasets by the parameter; and 

the merge and redundancy reducing function further comprises a selection function 
which applies a user-defined selection parameter whereby the subset is restricted relative to 
20; the sum of the datasets by the parameter, wherein the parameter is selected from the group 
consisting of source, author and pathway. 

In a third principal embodiment, the invention provides a computer-based system for 
creating a targeted collection of sequences from a dataset comprising sequence identifiers 
corresponding to natural complex biopolymer sequences and linked to corresponding first 
25 annotations, the system comprising: 

a) an integration function which merges the dataset with a database comprising second 
annotations attributable to and correlated with at least a subset of the sequence identifiers or 
sequences of the dataset and which links the second annotations to the corresponding 
sequence identifiers of the subset; and 
30 b) a tabulation function which creates and outputs the targeted collection of sequences 

in the form of a data table comprising, configurable by and sortable by the sequence 
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identifiers of the subset and the second annotations. 

The system may optionally incorporate the following limitation: 
the second annotations comprise data attributable to and correlated with at least a 
subset of the sequence identifiers or sequences of the dataset, said data selected from the 
5 group consisting of: gene expression data, sequencing data, genotype data, polymorphism 
data and clinical data. 

In yet another embodiment, the invention provides a computer-based system 
incorporating the elements of the first, second, and optionally, the third principal 
embodiments described herein. 

10 In a particular embodiment, the recited systems and methods have been implemented 

in a computer tool called ARROGANT. This program has been developed to facilitate the 
identification, analysis and comparison of collections of genes or clones. ARROGANT, in the 
analysis mode, is a comprehensive tool for providing annotation to large gene collections. 
ARROGANT takes in a large collection of gene identifiers and associates it with other 

15 information collected from many sources like sequence annotations, pathways, homology, 
polymorphisms, artifacts etc. to help the researcher draw scientific conclusions, 
understanding, and proceed with future experiments. The simultaneous annotation for a large 
assembly of genes makes the collection of genomic / EST sequences truly informative. For 
example, if the collection of genes is used for microarrays, ARROGANT predicts 

20 cross-hybridization with the members on the array and the entire UniGene database to help 
the researcher to design probes that avoid cross-hybridization or alerts the user of their 
presence. In the design mode, ARROGANT assists in compiling a gene collection, using 
several different databases simultaneously, queried with keywords and their synonyms. 
ARROGANT, in one integrated package, also facilitates the design of expression / 

25 resequencing microarrays by designing primers, looking for commercially available clones 
and designing probes for resequencing. The package also has a third mode of operation to 
eliminate sequence redundancies and duplicates from multiple gene collections. This is very 
useful in identifying redundancies due to sequences or clones having different accession 
numbers but representing fragments of the same gene. This simplifies comparing experiments 

30 from various research groups. ARROGANT has been successfully applied to many large 

gene collections for microarrays, complex multigenic trait projects, polymorphism discovery 



UTSD:668 



projects etc. 



BRIEF DESCRIPTION OF THE FIGURES 
Fig 1: Distribution of ARROGANT over three separate servers. 
Fig 2: Role of ARROGANT in designing expression/resequencing microarrays. 
Fig 3:Role of ARROGANT in analyzing microarray results. 
Fig 4: Flowchart for method for estimation of cross-hybridization. 
Fig 5: Block diagram of ARROGANT in design mode. 
Fig 6: Input page of ARROGANT in design mode. 
Fig 7: Sample output of ARROGANT from keyword search. 
Fig 8: Flowchart for primer design. 
Fig 9: Block diagram of ARROGANT in analysis mode. 
Fig 10: Selection of fields and selection of sort parameters page. 
Fig 1 1 : Final display page for analysis mode. 
Fig 12: Example output of Pathways. 
Fig 13: Schema for GenBank database 
Fig 14: Schema for UniGene database 
Fig 15: Schema for LocusLink database 
Fig 16: Schema for KEGG database 
Fig 17: Schema for HomoloGene database 
Fig 18: Schema for Research Genetics clone database 
Fig 19: Relationship across different databases. 
Fig 20: Implementation of the design mode. 
Fig 21 : Flowchart for the design mode. 
Fig 22: Implementation of the analysis mode. 
Fig 23: Implementation of the merging gene collection mode. 
Fig 24: Input for the merging gene collection mode. 
Fig 25: ARROGANT retrieves sequence properties 
Fig 26: ARROGANT sorts annotation table 
Fig 27: ARROGANT compiles unique list 
Fig 28: ARROGANT uses array elements to design probes 
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Fig 29: ARROGANT automates the selection of array elements 
Fig 30: ARROGANT eliminates sequence redundancies 
Fig 31a-b: ARROGANT tags false-positives 



5 DETAILED DESCRIPTION OF PARTICULAR EMBODIMENTS OF THE INVENTION 
The following descriptions of particular embodiments and examples are offered by 
way of illustration and not by way of limitation. 

The invention relates to computer-based systems and corresponding methods for the 
design and analysis of complex biopolymer sequence arrays. The recited biopolymer 

10 sequences are polypeptide or polynucleotide sequences, preferably natural sequences, and the 
arrays may be conceptual, digital arrays or tangible, molecular arrays (i.e. solid phase 
biochips). While exemplified with a server and web based application, the functionalities of 
the systems may be implemented in any convenient algorithm on any convenient platform, 
stored in any convenient medium, and accessed through any convenient port and interface. 

1 5 The functionalities of the recited systems may be independently segregated or two or more 
functionalities may be implemented together in a single process. Similarly, in the 
corresponding methods, the steps may be independently segregated or two or more steps may 
be effected contemporaneously. Of course, the recited systems and methods may further 
comprises further functionalities and steps to accomodate user preferences and database 

20 requirements and these may be implemented in concert or independently of those recited. 

In a first principal embodiment, a design mode operation, the invention provides a 
computer-based system for creating a targeted collection of sequences from a dataset 
comprising sequence identifiers corresponding to natural complex biopolymer sequences and 
linked to corresponding annotations. For example, in the design mode, the system can 

25 compile a collection of genes and/or proteins relevant to a user-determined study. The system 
is amenable to any selected nomenclature of sequence identifiers, such as accession numbers 
and a wide variety of corresponding annotations, which may include such attributes as source 
species, author, deposit date, associated diseases and pathways, polymorphisms, length, etc. 
Similarly, the system is amenable to any convenient dataset, unconstrained by size, 

30 complexity or public availability. Exemplary datasets include locally maintained datasets, 
GenBank, Medline, KEGG, etc. Furthermore, the system can implement a plurality of 
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datasets, preferably at least three, more preferably at least five, simultaneously or 
sequentially. 

In this first embodiment, the system effects four functionalities. The first 
functionality is a search function which searches the annotations of the dataset according to a 
5 user-defined criterion and outputs a first subset of the dataset restricted by the criterion. The 
user may select any criterion encompassed by annotations of the dataset(s), including 
keyword and more concept-based search criteria. Furthermore, the criterion is preferably one 
of a plurality of user-defined criteria, preferably at least three, more preferably at least five, 
and may include multiple keywords, such as "cardiac disease", and preferably uses logical 

10 operators like 'and/or", and the search function searches the annotations of the dataset 

according to all the criteria and outputs a first subset of the dataset restricted by the criteria. 

The second functionality is a redundancy reducing function which compares the first 
subset with a first database correlating the sequence identifiers of the first subset with 
syngeneic biopolymers and outputs a second subset of the dataset having reduced unique, 

15 natural complex biopolymer redundancy relative to the first subset, and preferably, which 

substantially eliminates unique, natural complex biopolymer redundancy relative to the first 
subset (i.e. wherein the second subset of the dataset has no unique, natural complex 
biopolymer redundancy). The database need not be perfectly syngeneically ordered (i.e. 
natural polynucleotides or polypeptides ordered by common source gene) and the recited 

20 correlation need not be perfect, but only sufficient to effect the required reduction in 
redundancy. Exemplary databases include locally maintained syngeneically ordered 
databases, UniGene and LocusLink. Preferably, the database is one of a plurality of 
databases, preferably at least three, more preferably at least five, correlating the sequence 
identifiers of the first subset with syngeneic biopolymers, and the redundancy reducing 

25 function compares the first subset with the databases and outputs the second subset of the 
dataset. 

The third functionality is a selection function which applies to the second subset a 
user-defined selection parameter and outputs a third subset restricted relative to the second 
subset by the parameter. As with the recited criterion, the used may select any parameter 
30 encompassed by annotations of the dataset(s), including source, species, author and pathway 
parameters. Furthermore, the paramter is preferably one of a plurality of user-defined 
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selection parameters, preferably at least three, more preferably at least five, and preferably 
uses logical operators like 'and/or", and the selection function applies to the second subset the 
parameters and outputs the third subset restricted relative to the second subset by the 
parameters. 

5 The fourth functionality is a tabulation function which creates and outputs the targeted 

collection of sequences in the form of a data table comprising, configurable by and sortable 
by the sequence identifiers of the third subset. The data table generally includes a variety of 
information related to each sequence identifier. The data table may be exported in any 
convenient format, preferably a format that is platform-independent, web-browser compatible 
10 and/or widely used, e.g. text, tab delineated, commercial database format, e.g. Excel, Access, 
etc., etc. 

This first principal embodiment preferably also comprises an expansion function 
which searches a second database for synonyms of the sequence identifiers of the first, second 
or third subset; note that the first database may be the same as the dataset. 

15- In a second principal embodiment, a comparison mode operation, the invention 

provides a computer-based system for creating a targeted collection of sequences from a 
plurality of datasets comprising sequence identifiers corresponding to natural complex 
biopolymer sequences, the system comprising two functionalities. For example, in the 
comparison mode, the system can resolve overlapping collection members, even if they have 

20 different identifiers for the same object. The system is amendable to comparing any of a wide 
variety of empirically, conceptually or computationally-derived databases and may operate to 
groupwise (e.g. pairwise) or simultaneously compare at least two, preferably at least three, 
more preferably at least five different datasets. 

The first functionality of this comparison mode of operation is a merge and 

25 redundancy reducing function which compares the datasets with a database correlating the 
sequence identifiers with syngeneic biopolymers and creates a subset of the sum of the 
datasets having reduced, and preferably substantially eliminated unique, natural complex 
biopolymer redundancy relative to the sum. As with the design mode, the database need not 
be perfectly syngeneically ordered (i.e. natural polynucleotides or polypeptides ordered by 

30 common source gene) and the recited correlation need not be perfect, but only sufficient to 
effect the required reduction in redundancy. Similarly, the merge and redundancy reducing 
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function may further comprise a selection function which, analogously with the design mode, 
applies a user-defined selection parameter whereby the subset is restricted relative to the sum 
of the datasets by the parameter. 

The second functionality of the comparison mode is a tabulation function which 
5 creates and outputs the targeted collection of sequences in the form of a data table 

comprising, configurable by and sortable by the sequence identifiers of the subset. As with 
the design mode, the data table may be exported in any convenient format, preferably a 
format that is platform-independent, web-browser compatible and/or widely used, e.g. Excel, 
Access, etc. 

10 In a third principal embodiment, an analysis mode operation, the invention provides a 

computer-based system for creating a targeted collection of sequences from a dataset 
comprising sequence identifiers corresponding to natural complex biopolymer sequences and 
linked to corresponding first annotations. As with the design mode, this system is amenable 
to any selected nomenclature of sequence identifiers, such as accession numbers and a wide 

15 variety of corresponding annotations, which may include such attributes as source species, 
author, deposit date, associated diseases and pathways, polymorphisms, length, etc. 
Similarly, the system is amenable to any convenient dataset, unconstrained by size, 
complexity or public availability. Furthermore, the system can implement a plurality of 
datasets, preferably at least three, more preferably at least five, simultaneously or 

20 sequentially. 

The first functionality of the analysis mode of operation is an integration function 
which merges the dataset with a database comprising second annotations attributable to and 
correlated with at least a subset of the sequence identifiers or sequences of the dataset and 
which links the second annotations to the corresponding sequence identifiers of the subset. 

25 The second annotations may be any of a wide variety of empirically, conceptually or 

computationally-derived annotations, such as data attributable to and correlated with at least a 
subset of the sequence identifiers or sequences of the dataset. Any measurable data may be 
used, including gene expression data, sequencing data, genotype data, polymorphism data and 
clinical data. 

30 The second functionality of the analysis mode of operation is a tabulation function 

which creates and outputs the targeted collection of sequences in the form of a data table 
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comprising, configurable by and sortable by the sequence identifiers of the subset and the 
second annotations. As with the design and merge modes, the data table may be exported in 
any convenient format, preferably a format that is platform-independent, web-browser 
compatible and/or widely used, e.g. text, tab delineated, commercial database format such as 
5 Excel, Access, etc., etc. 

In yet another embodiment, the invention provides a computer-based system 
incorporating the elements of the first, second, and optionally, the third principal embodiments 
described herein. In an alternative embodiment, the invention provides a system for the 
design and analysis of complex biopolymer sequence matrices, said system comprising: 
10 (a) a sequence input function which accepts user-defined sequences; 

(b) a sequence criteria input function which accepts user-defined instructions for 
sequence search criteria according to function or structure; 

(c) a sequence search function which operates on said criteria to search biopolymer 
; sequence databases according to said criteria to identify a population of diverse yet 

1 5 functionally or structurally restricted sequences; 

(d) a sequence population processing function which expands said population to 
encompass synonymous sequences and condenses said population by combining into common 
elements sequences which are syngeneic; 

(e) a matrix output function which presents results of said processing as a matrix of 
20-' elements corresponding to unique, syngeneic sequences; 

(f) a matrix data input function which accepts user-determined matrix sequence data; 

(g) matrix sort, search and filter input functions which accept user-defined instructions 
for sorting, searching and filtering matrix sequences; 

(h) a sequence matrix sort function which operates on said sorting instructions to sort 
25 matrix sequences; 

(i) a sequence matrix search function which operates on said searching instructions to 
search matrix sequences; 

(j) a sequence matrix filter function which operates on said filtering instructions to 
filter matrix sequences; 
30 (k) a report generating function which provides reports of the filtered matrix 

sequences; and 
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(1) a graphical user interface which provides fields for user inputting of the user- 
defined sequences, the search and processing criteria, the user-defined instructions for the sort, 
search and filter input functions, and for user visualization of the reports. 



5 EXAMPLE 

This design, comparison and analysis modes described above have been implemented 
in the computer tool known as ARROGANT. ARROGANT, in one integrated package, 
assists in compiling a gene collection starting from keywords and their synonyms, designs 
primers, looks for commercially available clones and designs probes for resequencing. Its 

10 focus is to enable the researcher to simultaneously work with a large collection of genes, their 
information and experimental data. Once a gene collection has been refined by the researcher 
(the design mode), ARROGANT then annotates the collection from many databases (the 
annotation mode). Experimental data can then be analyzed by association with this 
information (sequence annotations, pathways, homology, artifacts etc.) to aid in drawing 

15 scientific conclusions. In a large gene collection, sequence redundancies and naming 

differences can be resolved by ARROGANT (the collection comparison mode). ARROGANT 
calls upon the resources of several computers and many databases to design, annotate and 
compare collections. Once this computation is done, the collection can be manipulated, and 
experimental data can be integrated with it on a PC-based database server via a web server. 

20. (a) Features/Capabilities 

- Finds possible gene collection candidates by doing a keyword search on multiple databases. 

- Provides a comprehensive annotation for large list of genes. 

- The annotation is presented as a table, which can be downloaded and used for later reference; 
thus there is no need to follow any links to gather annotation. 

25 - The annotation table is searchable and sortable. 

- Experimental microarray, sequencing or other data can be overlayed on the annotation table. 

- Eliminates sequence redundancies and duplicates by eliminating sequences having different 
accession numbers but representing fragments of the same gene. 

- Combines several different lists of genes and provides the set of unique genes. 

30 - Designs primers for a large list of genes and lets the user design primers either in the 3' end 
or in the random region. 
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- Looks for commercially available clones. 

- Designs oligonucleotide probe sequences for resequencing microarrays. 

- Estimates the amount of cross-hybridization expected for each probe. 

(b) Hardware Requirements 

5 - Compute Server: Hardware-independent; run under HP/UX 1 1.00. 

- Database Server: Requires Windows NT/2000 and SQL Server 7.0. 

- Web Server: Win NT/2000 and IIS. 

(c) Software Requirements 

- Internet information web server 
10 - ASP 

- SQL Server 7.0 

- Visual basic 6.0 

- perl 5.005 or higher 

- BLAST (works with both WU-BLAST and NCBI-BLAST) - a code that identifies DNA or 
15 protein sequence similarity. 

- PRIMO (UTSW) - a code that designs DNA primers for laboratory experiments. 

- NCBI Tools (to get a fasta file for any accession number and to make a file BLASTable) 

(d) Databases Implemented 

- GenBank: DNA sequence database. 

20 - UniGene: database of clusters of DNA sequence that approximate those of genes. 

- LocusLink: database that links individual genes, their synonyms and all individual DNA 
sequence files that are associated with them. 

- KEGG Genome and Pathways: database of connections between networks of genes. 

- Research Genetics Clone Database: database of the commercially available clones. 

25 - Repetitive Element Database (Rep-X): database of potential polymorphic repeat sites that can 
be causative for disease. 

- HomoloGene: database of genes that have homologies in other species. 

Section 1 : Introduction to ARROGANT 
30 ARROGANT is a database driven tool developed to compile, annotate and merge large 

gene collections. NCBI, KEGG, Research Genetics and other custom databases have been 
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implemented locally since they were the most commonly used databases and were found to 
extensively cover various items of information related to sequences. The local implementation 
of various databases and tools (e.g. PRJMO, BLAST) makes ARROGANT independent of 
other applications and significantly improves its performance. The modular design facilitates 
5 addition of new databases with relative ease. ARROGANT has three modes of operation: 1 . 
Design mode (http://arrogant.swmed.edu/indexl.asp) 2. Analysis mode 
(http://arrogant.swmed.edu/index2.asp) 3. Merge gene collections mode 
(http://arrogant.swmed.edu/index3.asp). The design mode includes keyword searching for 
compiling gene collections and helps in the design of expression and/or resequencing 

10 microarrays. ARROGANT facilitates the design of resequencing and/or expression 

microarrays by looking for commercially available clones, designing primers and designing 
probes for resequencing. The analysis mode annotates large gene collections and estimates 
cross-hybridization for microarrays. When used for microarrays, ARROGANT takes over 
where ratios or clustering of sequences finishes to provide important data about genes and 

15 enables researchers to get a global view. ARROGANT has been used to pre-compute 

annotation for a large number of gene collections (http://arrogant.swmed.edu/precompute.asp), 
and the results are stored in the database. This allows quick retrieval of the data and lets the 
researcher dynamically sort the annotation table. The merging gene collection mode is used to 
avoid duplicates and redundancies in collections. ARROGANT provides a web based interface 

20 and hyperlinks various fields displayed in all the three modes. 



Section 2: ARROGANT Features 

1 . Compiles large gene collections starting from keywords: In a typical keyword search 
the researcher has to go to several different websites to look for genes of interest. After getting 

25 separate lists from different websites, she must combine them and avoid duplicates. 

ARROGANT facilitates this process by implementing several databases locally and providing 
simultaneous access to all of them. The databases may include major publicly available 
databases like GenBank, Unigene, polymorphism databases, clone databases, etc., and the 
keyword search may be done using Boolean operators like 'AND'/ 'OR'. 

30 2. Annotates large gene collections: One of the main features of the analysis mode is 

to gather all the relevant information from different databases for a large collection of genes 
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and present it as a searchable and sortable table rather than merely providing links to them. 
This table may be downloaded as a Microsoft Excel spreadsheet. The user is able to sort the 
data on different fields to enable her to select only a part of data to be viewed having certain 
characteristics, permitting biologically significant observations. For example, the data may be 
5 sorted based on repeats followed by map location to find some correlation between the ones 
having repeats on the array and their position in the genome. 

3. Designs primers for a large set of genes: Once a large collection of genes is 
compiled it is useful to design primers (to amplify sequences using PCR) for them 
simultaneously. The parameters for design of the primers may be set by the user, which may 

10 include melting temperature, oligo length and number of primers in each direction. 

Additionally there is an option to automatically design primers in 3 r end or random exon 
region. 

4. Finds commercially available clones: Typically spotted microarrays consists of 
either cDNA clones and/or PCR products. It is generally cheaper to buy the clones for a large 

15 collection of genes than to amplify them by PCR. ARROGANT is able to search for 

commercially available clones and the search for clones may be based upon UniGene cluster 
identifiers to avoid redundancies. 

5. Designs probes for resequencing: Resequencing microarrays consist of short 
sequences of 20-25 bases called probes. ARROGANT is able to output resequencing probes to 

20 be directly used by the Digital Optical Chemistry (DOC) project to make oligonucleotide 
microarrays. 

6. Tags cross-hybridization artifacts for expression microarray: Cross-hybridization is 
caused due to sequences sharing regions of high sequence identity that leads to false positives 
on the array. Another feature of the tool is to computationally estimate the amount of cross 

25 hybridization for each member on the array. 

7. Integrates experimental data: The tool is able to input an additional decimal number 
to be associated with each sequence. The number may correspond to expression data for 
microarrays, purity of repeats for polymorphism studies, priority level for a gene in a 
collection, etc. 

30 8. Provides Web based interface: The program is available for the users to be accessed 

over the web without the need to install any additional software. The web based user-friendly 
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interface isolates the researcher from the complexities of the code and makes it a more widely 
used tool. 

9. Adds hyperlinks: The data presented by the tool may be hyperlinked to its source. 
For example, accession numbers may be hyperlinked to their GenBank annotation. Although 
5 basic ARROGANT provides extensive coverage to various items of information, adding 
hyperlinks enables the researcher to view thier original source. 

Section 3: Materials and Methods 

3.1 Computational Tools: For optimal performance of ARROGANT, we distributed the 
10 project into three different components: ARROGANT required implementation of many 

different databases locally; the web-based interface communicates with the databases to store 
and retrieve data; ARROGANT implements other codes including PRIMO (to design primers), 
BLAST (to estimate cross hybridization) and FormatDB (to retrieve sequences for given 
GenBank identifiers). These three components (Database, Web-Interface and programs on 

15 HP/UX computer) were implemented on three dedicated computers for optimal performance 
as shown in Fig 1. The 'Database Server' has different databases implemented locally in a 
relational database format. 2. The 'Web Server' is used to host web pages for input and output 
of ARROGANT. 3. The programs running on HP/UX computer were utilized by developing 
an application in Microsoft Visual Basic 6.0 to communicate, run code (using telnet) and 

20 transfer files between the database server and the HP / UX computers. SQL Server 7.0 was 
chosen as the database software because of its ability to handle very large databases (having 
more than a million entries). Oracle was the only other database program capable of handling 
such large databases efficiently but is relatively expensive as compared to SQL Server 7.0. A 
web scripting application, Active Server Page (ASP) was used to develop various input and 

25 display pages. Perl programming language, especially useful for string operations, was used to 
develop other supporting scripts like parsing through BLAST output files to count hydrogen 
bond overlaps. The hardware and software configuration details are further detailed below. 

3.2 Materials and Methods for microarray experiment: For the utilization of ARROGANT in a 
30 microarray experiment, the steps involved are outlined in Figs 2 and 3. The first step in the 

design of a microarray experiment is to select members to be included on the array. The 
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researcher can search various databases like GenBank, dBEST, UniGene, etc. separately and 
combine into one list. ARROGANT helps in this process by enabling searching of candidate 
genes on several databases simultaneously. Once the members of the array are finalized, the 
next step is to look for the available clones or design primers for PCR. The output of a 
5 microarray experiment, genes correlated with expression level data, is usually fed into a 

clustering program. The researcher, based on her experience, has to interpret the results by 
individually going through different genes which are highly overexpressed or underexpressed. 
ARROGANT helps the researcher interpret the results by annotating all the genes in a tabular 
format and sorting the genes based on various items of information including expression level 
10 data. 

3.3 Materials and methods PGA (Program for Genomic Application): PGA is a nationwide 
attempt to elucidate the basic mechanisms and cellular responses underlying injury and 
inflammation and to identify functionally significant polymorphisms in human genes. 

15 ARROGANT may be used to annotate all the genes in the PGA list to point out all the 
polymorphic repeats and at the same time provide various important information (e.g. 
pathways, map location, etc.) about each gene in the list. ARROGANT incorporates 
additional candidate genes assembled from literature, collective local expertise of participating 
investigators and from other PGA centers. ARROGANT will also be used in the design mode 

20 to look for new candidate genes compiled by participating investigators. The new candidate 
genes suggested will again be annotated as one table in the analysis mode of ARROGANT. 
Another application identifies informative DNA sequence polymorphisms with respect to 
inflammation and remodeling processes within the heart and vasculature by using genetic 
association studies in large population groups, supplemented by family studies. 

25 

3.4 Materials and methods for cross-hybridization analysis: Yeast microarrays were used to 
study and quantify the effects of cross hybridization, as the entire yeast genome was available 
on a microarray. Each ORF (open reading frame) was BLASTed against every other ORF. The 
BLAST results were analyzed and the ORF's were ranked based on their sequence similarity to 

30 other ORF's. YBR1 12C was chosen as the target to be hybridized as it had the maximum 

sequence similarity to other ORF's. The results obtained from micro-array experiments were 



20 



UTSD:668 



then analyzed to find co-relation between the amount of cross hybridization and the BLAST 
output to determine the threshold for cross hybridization computationally. Various numerical 
scores were used to correlate the BLAST output with the expression level data obtained from 
the experiments. These were: 1. BLAST core, 2. Contiguous basepair overlap, and 3. 
Contiguous hydrogen bond overlap. We observed that cross hybridization becomes significant 
if the number of contiguous hydrogen bond overlaps was more than 65. This was used by 
ARROGANT as a threshold to identify potential cross hybridization. The algorithm used by 
ARROGANT for identifying cross hybridization is illustrated in Fig 4. Arrogant calculates 
cross hybridization in the analysis mode. A FASTA file for each gene on the array is obtained. 
Each sequence is BLASTed against every other sequence on the microarray and separately 
with the entire UniGene database. The results are tabulated as shown in Table 3.4.2 (two such 
tables are generated, one for potential cross hybridization by genes within the microarray and 
the other with the entire genome). 



Gene 


Number of 


Potentially cross- 


Maximum 


Sequence 




potentially cross- 


hybridizing genes 


overlap 


causing 




hybridizing genes 






maximum 
overlap 


<identifier> 


<#> 


<identitiers> 


<#bp> 


<sequence> 



Table 3.4.2: Output table for the cross-hybridization analysis in the analysis mode. 



3.5 Design Mode: The intent of the design mode is to let the user look for potential candidates 
associated with keywords to compile large gene collections and also help in the design of 
resequencing / expression microarrays . Fig 5 shows a block diagram of ARROGANT in 
design mode. The keyword search lets the researcher look for potential candidates in several 
different databases simultaneously. Fig 6 shows a snapshot of the input page for the design 
mode. This page serves to collect the following information from the researcher: 

1 . Email: Mandatory. Used later to email the results (primers, clones, FASTA files, 
probes etc.) to the researcher. 

2. Array Name: Mandatory. Used for the convenience of the user having more than one 

array. 

3. Keywords: List of Keywords to be searched. Logical operators like AND' / 'OR' can 
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be used. 

4. Select Databases: (At least one must be selected) Multiple databases maybe selected 
at one time; options include: a. GenBank; b. UniGene; c. LocusLink; d. KEGG; e. Research 
Genetics clone database (http://www.resgen.com). Any combination of the above databases 

5 may be used. 

5. Include File: (Optional ). The researcher can add a list of accession number tabs 
delimited with a floating-point number (any number in general which may indicate priority 
level, purity of repeat, expression level data, etc.). The 'include file' accession numbers are 
included in the final list irrespective of whether the keyword search finds it. These accession 

10 numbers are added to the list selected by keyword search before proceeding with the next step 
of designing primers / looking for clones or designing resequencing probes. 

3.5.2 Keyword Search: The search is done in such a way that the entries common to different 
databases are represented only once. Multiple keywords can be used in the search using 

15 Boolean operators like AND, OR. e.g. cancer AND aging, klotho OR kl. If no operator is 
specified then 'and' is assumed, e.g. 'aging klotho' is translated as 'aging and klotho'. The 
search is case-insensitive. The results are separated into pages and the user can go through 
different pages using the "Next" and "Back" button or clicking directly on the hyperlink of the 
page number. The user has the option to select a few entries from the keyword search or enter 

20 new keywords and move between pages until he finalizes the gene collection. Fig 7 shows the 
sample output obtained by doing a keyword search. The number of entries on each page can be 
entered by the user. 

3.5.3 Fields Retrieved: In general, a unique identifier, a short description and organism are 
25 retrieved from each database. Accession number is retrieved as the unique identifier from the 

research genetics database and all the NCBI databases which include GenBank, UniGene and 
LocusLink. GenBank numeric identifier (NID) from the KEGG database and pathway number 
from the KEGG pathway database are retrieved as unique identifiers. The unique identifier 
from all the databases are hyper-linked to their respective annotation. The output of a keyword 
30 search is similar to that obtained on the NCBI website. 
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3.5.4 Resequencing: Resequencing is done to study the mutations and sequence variation in a 
DNA sequence. ARROGANT helps in the design of resequencing microarrays by calculating 
oligonucleotide probes. The Digital Optical Chemistry (DOC) technology may be used to 
build custom resequencing arrays. The DOC technology uses Digital Light Processors (DLP) 

5 chips (Texas Instruments) to create 'digital masks' dynamically. These 'digital masks' are 
created by directing the light reflected off the DLP. This is possible since every individual 
mirror on the DLP can be controlled independently. Consequently, the DOC system has the 
advantage of being able to manufacture custom microarrays for resequencing/ expression and 
the number of array elements can reach 2,000,000. ARROGANT creates a series of probes for 
1 0 each sequence to be put down on the chip for resequencing. The program generates 1 6 probes 
of 25 mer after taking different parameters like melting temperature into consideration. This 
can be directly used by DOC to make the resequencing arrays. 

3.5.5 Expression: Expression microarrays are primarily used to study the expression of 

1 5 thousands of genes simultaneously. This technique is very useful if used in different cycles of 
cell development. Comparison studies involve studying the expression of thousands of genes 
from normal and pathological subjects. Different color dyes are used to differentiate between 
the normal and pathological samples.Genes behaving differently are usually studied further. In 
the case of expression studies, typically clones or PCR products representing a particular gene 

20 = sequence are spotted down. Our lab also has the spotted microarray technology to develop 
expression microarrays. The DNA which is spotted is typically around 0.4 kilobases to 2 
kilobases. In the case of expression mode, ARROGANT lets the researcher either look for 
commercially available clones or design primers to create PCR products. 

25 3.5.6 Design of Primers: ARROGANT uses a code called PRIMO available at 

http://atlas.swmed.edu. The code has been successfully used to design primers for a large 
number of PCR reactions. PRIMO uses a 'regions' file to design primers and amplify the 
specified region of interest. ARROGANT creates the 'regions' based on the user's selection to 
design primers either in the 3' or random region. ARROGANT lets the user modify the 

30 parameters used for the design of primers. The parameters include: 1. Oligo length: Length of 
the primer to be designed which is typically around 20 bases. 2. Tm: Melting temperature to 
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be used for PCR reactions. 3. Number of primers to select (per direction) : Number of forward 
and reverse primers to select (default =1). Fig 8 shows a flowchart for primer design. 

3.5.7 Commercially available clones: A database of clones available with Research Genetics 
5 has been implemented. The UniGene cluster identifier is first obtained for each sequence and 
this identifier (instead of GenBank accession/nid ) is used to search for clones. The advantage 
of using UniGene identifiers is that it represents the cluster and all accession numbers 
belonging to it. This avoids redundancies. This approach eliminates the possibility of getting 
different clones for the same gene represented by separate accession numbers. Research 
10 Genetics provides its own software called CMiner to find available clones. Results of CMiner 
have been compared to ARROGANT and found to be the same. 

3.6 Analysis Mode: The intent of the analysis mode is to annotate a large gene collection and 
display the annotation table in a format conveneint for the user to view. 

15 

3.6.1 Block Diagram: ARROGANT accepts a list of accession number tabs delimited by a 
floating point number. Using the various databases implemented locally it annotates the list of 
accession numbers with 36 different items of information. Fig 9 shows the block diagram of 
ARROGANT in the analysis mode. The researcher then has a choice of selecting the fields he 

20 would like to view. The researcher can also choose the fields on which the output table is to be 
sorted. The researcher can sort and sub-sort based on five different fields at one time. Fig 10 
shows the 36 different items of information to choose from and also the five different choices 
the user has to sort and sub-sort the table. 

The input page consists of the following fields: 

25 1. Email: Mandatory. The link to the results page (generated after the analysis mode is 

completed) is sent to the user. 

2. Array Name: Mandatory. Used for the convenience of the user having more than one 

array. 

3. Include File: Mandatory. The user can add a list of accession numbers tab delimited 
30 with a floating-point number to input a list of accession numbers to be analyzed. 

A count of number of entries found for each field is provided where the data from the 
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first sort column becomes "Not found / Unknown". An additional summary is also provided at 
the end of the table. For example, consider that the annotation table was first sorted based on 
'Homolog' field. At the point when the homologs are not found there would be a summary of a 
count of entries for each field and also at the end of the annotation table. So if the user would 
5 want to do a mouse 'knock-out' experiment he would have a summary of a count of each field 
for only the ones having homologs. Fig 1 1 gives the snapshot of the output display page in the 
analysis mode. 

3.6.2 Fields Retrieved: ARROGANT in the analysis mode annotates each entry in a gene 
10 collection with 36 different items of information. The following section explains the 
significance of each field and its source of retrieval. 

1 . Accession: The accession number is the most widely used identifier to represent 
sequences. Accession numbers do not change even if information in the record is changed at 
the author's request. Consequently, accession numbers are used to input the list of gene 

15 collection to ARROGANT. Examples: NT_123456 (constructed genomic contigs); 
NM_123456 (mRNAs); NP_123456 (proteins); NC_123456 (chromosomes). 

2. Unique Identifier (NID) : "Genlnfo Identifier" (GI) is another sequence 
identification number used by GenBank . If a sequence changes in any way, a new GI number 
is assigned. ARROGANT retrieves NTD for each GenBank accession number, e.g. 4557722. 

20- 3. Related Sequences: All sequences belonging to the UniGene cluster of the current 

entry are considered related and retrieved. Related sequences are obtained from the NCBI 
UniGene database. 

4. Taxonomy: Taxonomy of the organism to which the given sequence belongs is 
obtained from NCBI GenBank database implemented locally. 
25 5. Repeat: Potentially polymorphic repeats for the given sequence cluster are identified. 

This is retrieved from the output generated by our code Rep X. 

6. Hairpin/Palindrome: Hairpins and Palindromes are found in the sequence by Rep X. 

7. Homology: HomoloGene (Zhang et al, J. Comp. Biol. 2000) database is used to find 
homologs/ orthologs. UniGene identifiers are used to represent the whole cluster. 

30 HomoloGene database calculates homologs by nucleotide sequence comparison between all 
UniGene clusters for each pair of organisms. Homology information might be useful for 
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'knock-out' experiments. 

8. Research Genetics Clone ID: The researcher might want to know wheather a clone is 
commercially available for a given sequence. ARROGANT outputs clone identifiers available 
with Research Genetics. The clone query is done using its corresponding UniGene cluster 

5 identifier. 

9. IMAGE: ARROGANT also lets the researcher retrieve information about whether a 
clone can be obtained from the IMAGE consortium. 

10. Cross Hybridization: This field is specific to using ARROGANT for expression 
profiling microarrays. An artifact sometimes observed in the results obtained from an 

10 expression profiling microarray experiment is that some sequences might hybridize to other 
sequences to which they are significantly similar. ARROGANT points out the potential 
cross-hybridizing candidates due to regions of sequence similarities. It tags the potential 
candidates and lists a score of the total false signal expected and the corresponding 
cross-hybridizing genes. This information can also be used by the researcher to design primers 

15 so as to exclude the regions causing cross-hybridization. 

1 1. cDNA Source: Determines the source of DNA for the particular sequence. It is 
obtained from the "cDNA Source" tag of UniGene. This helps the researcher know the 
possible places where the gene was found to be expressed. 

12. GDB ID: The Genome database is an important resource used to retrieve 

20 information about a sequence. ARROGANT retrieves the GDB identifier and hyperlinks the 
identifier to its original source. This is obtained from LocusLink database. 

13. Gene Name: A given sequence can have a number of different gene names. 
ARROGANT retrieves the official gene names for a particular sequence or of the UniGene 
cluster to which it belongs. Retrieved from UniGene as the 'Gene' tag, LocusLink as the 

25 'Official Gene Name', Genome Database as 'Description' tag and the KEGG database as the 
'Name' tag. 

14. Expression Data: This is an additional field entered with the accession numbers, 
which could represent expression level data / log expression value (for microarray data), purity 
of repeat (polymorphism studies), priority of genes on the list or order of genes on the list. The 

30 expression level data can be a number anywhere from 0.00001 to 1000000. 

15. Gene function: ARROGANT provides a short summary of the function or 
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description of the sequence / cluster of the current sequence. This is retrieved from UniGene as 
the 'Title' tag and from LocusLink as the 'Summary' tag. This helps the researcher to get an 
idea of the possible function for a given sequence. 

16. Synonyms: A list of symbols and names used for a given gene sequence. It is 
5 retrieved from LocusLink as 'Official Symbol' tag. This can be used by the researcher to 

identify all the names/symbols associated with the given sequence. These could be further used 
for a keyword search. 

17. Pathways: Pathways that the gene sequence might be involved in. It is retrieved 
from the KEGG database. Fig 12 shows an example of a pathway displayed by ARROGANT. 

10 The pathway information can help the researcher deduce important inferences. For example in 
a microarray studies, the expression level data when correlated with pathways can help include 
new genes in a pathway or associate new pathways with diseases and phenotypes. 

18. SNP: Single Nucleotide Polymorphism records for a particular sequence are 
retrieved from the LocusLink database. This might help associate mutations with specific 

15 pathological conditions. 

19. Title: Brief description of sequence obtained from GenBank database which 
includes information such as source organism, gene name/protein name, or some description 
of the sequence's function (if the sequence is non-coding). 

20. PubMed: Provides a link to the various articles and journals related to the sequence 
20 " from the PubMed database. PubMed entries are retrieved from LocusLink database. This helps 

the researcher view an abstract of the article describing the function of the particular gene 
sequence. 

21. Reviewed RefSeq: ARROGANT outputs reference sequences based on their 
LocusLink cluster. This helps eliminate different accession numbers representing the same 

25 gene. 

22. UniGene ID: Represents the identifier of the UniGene cluster to which the given 
sequence belongs. The UniGene identifier is further linked to its UniGene annotation which 
can be a useful source of information to the researcher. 

23. LocusLink ID: Represents the identifier of the LocusLink cluster to which the 
30 given sequence belongs. This is again hyperlinked to its LocusLink annotation. Many 

researchers might frequently want to convert their lists from one identifier to another. Using 
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the annotation mode of ARROGANT, accession numbers can be converted to LocusLink, 
UniGene, and KEGG identifiers. 

24. Alias Symbol: ARROGANT lists alias symbols used for a given gene sequence. 
They are retrieved from LocusLink using Alias Symbol tag. 
5 25. Alias Protein: ARROGANT lists alias proteins for a given gene sequence. They are 

retrieved from LocusLink using 'Alias Protein' tag. 

26. Phenotype: Phenotypes found and associated with mutations in this gene. They are 
obtained from the LocusLink tag 'Phenotype 1 . 

27. Phenotype ID / OMPM ID: For humans, this represents OMLM (Mendelian 
10 Inheritance in Man) number. They are retrieved from LocusLink as 'Phenotype ID' tag. 

28. Map Location: It indicates the location of the gene on the chromosome map. Map 
location is retrieved from the 'Map' tag of LocusLink and also from its UniGene entry. 

29. Map Link: Provides a link to the graphical representation of the gene on the 
chromosome. This again is retrieved from LocusLink- 'Map' tag. 

15 30. Map Type: Map-Type is the type of map information (G = genetic, 

C=Cytogenetics). This also is retrieved from LocusLink- 'Map' tag. 

31. STS - Markername: STS are sequence tag sites, which are short sequences that are 
operationally unique in the genome, used to generate mapping reagents. It is retrieved from the 
STS -Markername tag of LocusLink. 
20 32. STS - Chromosome: Retrieved from the STS-Markername tag of LocusLink 

33. STS - ID: Retrieved from the STS-Markername tag of LocusLink 
Each item is hyperlinked to the help file which exactly defines what each field means 
and how it is retrieved. 

25 3.6.3. Sorting: ARROGANT lets the user sort the final output display based on five different 
fields. The intent is to allow the user to look at the output in a convenient manner. This would 
help the user to answer questions like 'Which of the genes in this large gene collection have 
known SNP entries and also have a mouse homolog which could be used for a knock-out 
experiment?' Consider an example where the researcher has a large gene collection to study 

30 polymorphism. The researcher might first group the genes based on their chromosome location 
by sorting them based on chromosomes. Assume that the user is interested in looking for 
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polymorphic genes on chromosome 3. He could select the second level of sort as repeats 
which would bring the genes belonging to chromosome 3 and having polymorphic repeats on 
top of the list. Further the researcher might want to know for how many of these top genes are 
the pathways known. This can be achieved by selecting the third level of sort as pathways. 
5 Thus the user can sort and sub-sort to analyze the large gene collection in a convenient 
manner. 

3.6.4. Add Accession Numbers: ARROGANT lets the user add accession numbers to her 
already submitted list in the analysis mode. The user has the option of adding more entries to 
10 her list while viewing the annotation of her previous list. The new list thus created 

automatically avoids duplicates and sequence redundancies. This lets the user add more entries 
to her list assembled earlier. For example, the user could have annotated a list of four thousand 
genes associated with cancer using the analysis mode of ARROGANT and can add thousands 
more entries to be annotated and combined with this list years later. 

15 

3.7 Merging Gene Collection Mode: In merging gene collection mode, the user could submit 
two or more lists of genes which could be combined into one unique list avoiding all 
duplicates. The redundancy is avoided by eliminating accession numbers belonging to the 
same UniGene cluster. The input page to the merging gene collection mode consists of the 
2G_ following input fields. 1. User Email: Mandatory. The merged gene list and the original gene 
collection is sent to the user on this e-mail address. 2. Array name: Mandatory. The name to be 
associated with the gene collection. 3. Input file: Mandatory. Input one file containing 
sequence redundancies. 

25 Section 4: Implementation 

This section deals with the details of the software implementation of ARROGANT. 
ARROGANT is a database driven software. The section talks about databases implemented by 
ARROGANT, algorithms used in the design, analysis and the merging gene collection mode. 

30 4.1 Databases: As described, ARROGANT combines results from several different databases. 
All the databases are implemented in the relational database format in SQL server 7.0. In the 
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following section each database is introduced and its implementation is covered in brief. 



4.1.1 GenBank: GenBank, an annotated collection of all publicly available DNA sequences 
provided by NTH, is the biggest and the most used publicly available database (Nucleic Acids 

5 Research 2000 Jan l;28(l):15-8). There are approximately 10,897,000 sequence records as of 
February 2001 (http://ncbi.nlm.nih.gov). The complete release notes for the current version of 
GenBank are available at ftp://ncbi.nlm.nih.gov/genbank/gbrel.txt. The GenBank database is 
the single most important database to search for possible gene candidates. Each GenBank entry 
has a unique identifier called accession number. ARROGANT uses accession number as its 

1 0 primary key to link different databases. ARROGANT uses GenBank database in design and 
analysis mode. GenBank is implemented as a separate database on the server called 
'genelfullgenebank' to improve the speed performance, as the database is very large containing 
approximately 10.8 million entries. The database is implemented as a single table, see Fig 13. 
ARROGANT GenBank database implemented in SQL Server 7.0 does not include the actual 

15 sequence for each entry. This is obtained using the NCBI tools implemented locally on our 

HP/UX computers. A shell script 'getgb' compares files present locally with its original source 
on the web and downloads only the ones not existing or having a different file size from 
ftp://ncbi.nlm.nih.gov/genbank. The files are unzipped, combined into one huge file, split into 
smaller files of approximately equal sizes and then reformatted and can then be directly 

20 imported into the database using the 'bulk insert' script. 

4.1.2 UniGene: UniGene partitions GenBank EST sequences into a non-redundant set of 
gene-oriented clusters. Each UniGene cluster contains sequences that presumably represent a 
unique gene, as well as related information such as the tissue types in which the gene has been 

25 expressed and map location. The UniGene database was chosen to be a part of ARROGANT 
(see Fig 14) for the following reasons: 1. Avoid Redundancy: ARROGANT uses UniGene 
database to avoid redundancies by not including sequences having different accession numbers 
but representing the same UniGene cluster. ARROGANT uses this in the merging gene 
collection mode to combine different lists into one unique collection. 2. The UniGene database 

30 includes gene sequences as well as hundreds of thousands of expressed sequence tag (EST) 
sequences. 3. Additional Annotation: Provides additional annotation for a given gene 
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sequence, e.g. cDNA source, which is used to look for keywords (design mode) and annotate 
gene collection (analysis mode). As a result UniGene database is used in all the three modes 
by ARROGANT. Perl scripts combine similar files 

(ftp://ncbi.nlm.nih.gov/repository/UniGene/) of different organisms together, convert the files 
5 into various files of specific format which can be imported directly into the database tables 
using the import function in SQL Server 7.0. 

4.1.3 LocusLink: LocusLink is NCBI's attempt to integrate and provide a single query 
interface to clustered sequences and make available descriptive information about genetic loci. 

10 However, LocusLink does not provide annotation to a collection of genes. ARROGANT 

extends its capabilities by incorporating LocusLink database. Sequence accessions include a 
subset of GenBank accessions for a locus, as well as a new type, the NCBI Reference 
Sequence (RefSeq). LocusLink provides a reference sequence for each locus cluster. 
LocusLink database is used by ARROGANT in the design and analysis mode, see Fig 15. 

15 Series of Visual Basic executables import files into the database, downloaded from NCBI 
(ftp://ncbi.nlm.nih.gov/refseq/LocusLink/LL_tmpl). 

4.1.4 KEGG Genome and Pathway Database: ARROGANT not only combines different 
databases from NCBI but also uses the KEGG databases. Kyoto Encyclopedia of Genes and 

20 Genomes (KEGG) makes available, information pathways consisting of interacting molecules 
or genes by using the current knowledge of molecular and cellular biology (Kanehisa, M., 
Oxford University Press 2000). In addition KEGG database also provides additional 
annotation used by ARROGANT to look for keywords and annotate gene sequences. As a 
result KEGG database is used by ARROGANT in both design and analysis mode, see Fig 16. 

25 The files downloaded from KEGG (ftp://kegg.genome.ad.jp/genomes/genes/ ) are combined as 
one, split into smaller files and the Visual Basic executable is used to update the tables. A file 
containing additional pathway information is used 
(ftp://kegg.genome.ad.jp/pathways/map_title.tab). 

30 4. 1 .5 HomoloGene: The HomoloGene database provides homologs /orthologs, which is used 
as a field in the annotation of large gene collection by the analysis mode, see Fig 17. It 
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primarily uses the UniGene cluster identifier to search for homologs / orthologs. Accession 
numbers and LocusLink identifiers may also be used. HomoloGene uses nucleotide sequence 
comparisons to calculate orthologs and homologs, between all UniGene clusters by each pair 
of organisms. The HomoloGene database is downloaded from 
5 ftp://ftp.ncbi.nlm.nih.gov/pub/HomoloGene/hmlg.ftp. Perl scripts format the downloaded file, 
which is further imported into the database. Special character ,A ' is used as the delimiter to 
import the file into the database. 

4.1.6 Research Genetics Clone Database: Research Genetics commercially distributes selected 
1 0 clones from the IMAGE consortium. The catalog of clones available at Research Genetics can 

be downloaded at ftp://ftp.resgen.com/pub/sv_libraries/RG_Hs_seq_ver_101100.txt. The 
catalog contains annotation related to the clones like accession number, gene name, cluster ID, 
insert size, markers, etc. ARROGANT stores this catalog locally in the database, which is used 
to find commercially available clones and search for candidate genes in the design mode, see 
ISM Fig 18. 

4.1 .7 Relationship of different databases to one another: ARROGANT facilitates 
simultaneous keyword searches and annotates gene collections using different databases. This 
is achieved by implementing all the databases locally in SQL server 7.0 database. As shown in 

20 Fig 19 GenBank database is the only database linked to all other databases. GenBank is linked 
to KEGG database using NID (Numeric Identifier) and to all other databases by accession 
number. UniGene and LocusLink databases both contain LocusLink Identifier as well as 
UniGene identifier. The HomoloGene database is linked to GenBank by accession number, to 
UniGene by UniGene identifier and to LocusLink by LocusLink identifier. RepX database 

25 contains UniGene identifier. The KEGG database is also linked to LocusLink via LocusLink 
identifier. Thus ARROGANT integrates different databases from both NCBI and KEGG. 

4.2 Implementation of the Design Mode: ARROGANT in the design mode can be used to find 
candidate genes using keyword search, and design expression/resequencing microarray. The 
30 keyword search can be done online and the selected entries are queued in the database. Fig 20 
shows the implementation of ARROGANT in design mode. Keywords can be entered and 
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databases can be selected in the input page e.g. 'Cancer AND chromosome 3' can be searched 
in GenBank and LocusLink. The user can select a few genes of interest by going through 
different pages from the keyword search. Each gene is hyperlinked to its annotation. The user 
can also input new keywords and select from these new entries to add to the final list. The 
5 final gene list is a set of accession numbers/ GenBank unique identifiers. The final list is then 
inserted in a database table and the request is queued. The requests are processed one by one. 
Fig 21 shows the flow of events after the request is queued. The list of accession numbers is 
first transferred to the HP / UX computer where duplicates are eliminated. The FASTA files 
are obtained for these accession numbers using the NCBI tools (implemented locally), which 
10 run on all the NCBI sequence databases to make sure that none of the entries are lost e.g. 
Human EST sequences, patent sequences, non-redundant sequence database, etc. 

After the gene collection is finalized, the user can further use ARROGANT to design 
expression or resequencing microarray. For expression microarrays the user gets to choose to 

1 5 either design primers and/or find commercially available clones. If the option includes 

designing primers then depending on the whether the user wants it designed in the 3' / random 
exon region an appropriate regions file is created. Once the 'regions' file is created PRTJVIO is 
run on the combined FASTA file to design primers for each sequence so as to amplify the 
specified region. The output files are transferred back to the windows terminal. Depending 

20 upon the user selection, clone query may run on the list to find out the possible clones. Finally 
the request is cleared from the database and the user is sent the following attachments: 1. List 
of GenBank identifiers for genes selected. 2. Combined FASTA file for all the genes. For 
expression mode, 3. File containing designed primers for all the sequences. 4. A log file 
specifying the parameters used for the design of primers. 5. A file containing Research 

25 Genetics clone identifiers for clones available with Research Genetics. For resequencing 

mode, 6. A file containing resequencing probes for all the sequences on the oligonucleotide 
microarray. This summarizes the algorithm of ARROGANT in design mode. Consequently, 
ARROGANT in the design mode can be used to compile a large gene collection. For 
microarrays the design can be further extended to help generate expression / resequencing 

30 microarrays. 
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4.3 Implementation of the Analysis mode: ARROGANT in the analysis mode, is intended to 
be an annotation tool for large gene collections. The following section describes the user 
interface and the algorithm used in the analysis mode. Fig 22 summarizes the working of 
ARROGANT in the analysis mode. The list of accession numbers is first stored into the 
5 database table and the request is queued. The queue is processed one by one. The processing 
begins by running a stored procedure called 'testcode' which pulls the different fields of 
annotation related to the sequence as listed in the earlier section of block diagram description. 
ARROGANT then proceeds to estimate the amount of cross-hybridization possible within the 
chip and with the entire UniGene database as described in section 3.3. The results related to all 
1 0 the fields on the array are stored in a database and a link unique to the current array is 

generated using its taskid. This link is emailed to the user. Thus the user could click on the 
link and access her collection of genes. The user can further add more accession numbers to 
her collection by using the 'Add Accession' option on the final display page. 

15 4.4 Merging gene collection: The intent of this mode is to combine genes from different 

collections/microarrays into one unique list. The duplicates are avoided by first converting the 
accession numbers to their corresponding UniGene cluster identifiers and then retaining the 
ones which are unique. Fig 23 shows the algorithm used to merge different gene collections. 
Fig 24 shows the steps need to be taken. The user needs to first combine all the gene 

20 collections in one file. This is the input to the merging gene collection mode and the output 
contains only the unique ones. 

4.5 Tools for maintaining the databases: Tools have been developed to maintain the various 
databases used by ARROGANT. The tools include various Visual Basic executables and perl 
25 scripts which either directly update the database or create a file which can be directly imported 
into the database. The update of each individual database is described in its implementation 
section. The databases are maintained in a semi-automated fashion. The person in-charge of 
maintaining the program has to run scripts and later inspect the lock files to verify that the 
databases were updated properly. This is usually done every fifteen days. 

30 

Section 5: Applications 
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ARROGANT provides a useful tool for working with a large collection of genes and for 
design and analysis of microarrays. In addition, ARROGANT may be used in a wide variety 
of other applications: 

5 5.1 Retrieve properties to annotate large gene collections: Several different databases 

implemented locally enable ARROGANT to annotate any given sequence. In the analysis 
mode, ARROGANT accepts a list of accession numbers from the user and annotates it with 36 
different items of information. This is the only tool which annotates a large number of genes 
simultaneously and presents the data in a tabular format. The table can be downloaded locally 

10 and stored on the computer without having to follow any hyperlinks. Fig 25 summarizes the 

role of ARROGANT to retrieve properties for a large gene collection. For example, consider a 
researcher needing to know the chromosome location and the UniGene cluster identifiers for a 
collection of genes. As long as the list consists of only a small number of genes (less than 30 
- 40 ) she could do it manually by searching for individual genes on the web. With the recent 

1 5 advancements in technology and the completion of the Human Genome Project researchers are 
now dealing with hundreds and thousands of genes. ARROGANT offers unprecedented value 
for retrieving such gene properties for large gene collections. 

ARROGANT provides a tool of choice for many users wanting to annotate their large 
gene collection and have the information stored locally on their computer. This is the only tool 

20- the researcher needs to use in order to annotate completely a large list of gene collections. This 
saves the researcher from exploring multiple different websites and trying to compile all the 
information, which becomes practically impossible as the number of genes in the list 
increases. The details of each field used for annotation are described above. 

25 5.2 Discovery Tool: ARROGANT, when used in the analysis mode makes it easier for the user 
to generate biological inferences and recommendations for future experiments by enabling the 
user to view a large number of items of information, sorted based on multiple different fields. 
When used for microarrays, the analysis mode of ARROGANT associates the given gene 
sequence with 36 different items of information. The user can also overlay a floating point 

30 number associated with each sequence which may be expression level data (microarrays), 

priority level, purity of repeats (for polymorphism studies ), etc. The use of ARROGANT in 
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assisting discoveries can be illustrated with an example. Consider the user has entered her 
gene list with the expression level data in the analysis mode. After selecting the fields to be 
viewed, there are five levels of sorting available to the user. The user can view the genes 
expressed highly on top of the list by first sorting on expression level data. Further to get an 
5 indication of the genes expressed highly and belonging to a common pathway the second level 
of sort could be chosen as pathways. So if out of the ten highly expressed genes eight belong 
to a common/ related pathway, and pathways for the remaining two were unknown, then it 
might indicate that the remaining two belong to the same or related pathways. In this way 
ARROGANT can help answer questions like 'How many genes that expressed more then 5 
10 fold and having polymorphic repeats also have a mouse homolog which could be used for 

further knockout experiments?' which can be achieved by sorting on expression data, repeats 
and homologs. This is easily possible because the various items of information are in one place 
and the data can be sorted in various ways to filter the interesting ones on top of the list. Fig 26 
illustrates the use of ARROGANT to assist discoveries. 

15 

5.3 Compile a gene collection: keyword search, avoiding redundancies and duplicates. The 
merging gene collection mode can be used in combination with the keyword search to compile 

-= a unique collection of genes. Fig 27 shows how ARROGANT helps the researcher in 

compiling a gene collection. ARROGANT permits surfing through several different databases 

20 at one time. A parent schema of all these different databases has been implemented which lets 
the researcher select the various databases in any combination for a keyword search. This 
eliminates the need for the user to compile several different lists obtained from many different 
databases and trying to compile one unique list of genes. ARROGANT also avoids sequence 
redundancies by allowing only one gene per UniGene cluster. The keyword search can be done 

25 using logical operators like AND / OR', e.g. 'klotho or kl', 'aging and telomerase'. With the 
merging gene collection mode, the user can also chose to include an additional list of 
accession numbers of her genes (in addition to that obtained from the keyword search) to be 
included irrespective of whether the keyword search finds it, and a unique list is compiled 
automatically to avoid duplicates. 

30 

5.4 Design of oligonucleotide probes for resequencing: ARROGANT provides the user the 



36 



UTSD:668 



unique facility of designing probes after selecting the list of genes using a keyword search 
and/or adding genes from her own list. The probe design software has been tested and 
implemented for more than two years in our laboratory. The probes designed have been used 
to develop chips using DOC (Digital Optical Chemistry) technology developed in our 
5 laboratory. For each gene 16 probes of 25 mers each are designed taking parameters like 

melting temperature, oligo length etc. into consideration. Fig 28 illustrates how ARROGANT 
can be used to select members and to design oligonucleotide probes for resequencing arrays. 
Thus ARROGANT serves as a handy tool for design of resequencing arrays. The design 
includes selecting unique elements on the array (relevant to the topic of study) and designing 
1 0 probes to represent the sequences on the array. 

5.5 Design of expression microarrays: Expression microarrays often have elements of around 
0.4Kb- 2Kb length spotted; typically, they are either PCR products or clones commercially 
available. As shown in Fig 29, ARROGANT automates the selection of elements to be spotted 
15 on the array by first looking for possible elements by doing a keyword search and then for the 
ones selected, the user can either chose to design primers or to look for commercially available 
clones or both. 

5.5.1 Design of Primers: ARROGANT also automates the process of design of primers after 
20 selection of the sequences. The design of the primers can be made either in the 3' region or 

random exon region. Depending on the user selection, the 'regions file' which specifies the 
region for the design of primers is done so as to select the 3' region or random region. PRTMO 
(supra) designs primers based on the regions file and the conditions mentioned in the criteria 
file. The conditions for the design of primers can be modified by changing the criteria file. 
25 Oligo length, melting temperature for PCR and number of primers in each direction can be 
edited by the researcher. 

5.5.2 Commercially available clones: ARROGANT lets the researcher look for commercially 
available clones from Research Genetics, Inc. (Huntsville, AL). This can be done either by 

30 submitting the list of genes in the analysis mode and then selecting the Research Genetics 
clone column from the list of annotation. The researcher can also use the design mode of 
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ARROGANT to look for clones by choosing the expression microarray design. The results 
obtained have been compared to the output of CMiner (Research Genetics, Inc.) and results 
have been found to tally. 

5 5.6 Merging gene collections: In the merging gene collection mode the user has the 

opportunity of combining several different expert list of genes and obtain one unique list. For 
example consider that the researcher has three separate gene collections related to studying 
cancer. The first list might have been obtained by doing a keyword search on a popular 
database like GenBank. The other list might be obtained from the Cancer Genome Anatomy 
10 Project (CGAP, Cancer J 7(1)52-60, Schaefer et al, 2001) recommendations and the third 

maybe the genes she feels play a role in cancer. These three lists if combined manually by the 
researcher might contain a lot of different gene identifiers representing the same gene. She 
then manually needs to eliminate such duplicates by searching for each gene. ARROGANT 
automates this process by allowing only one gene per UniGene cluster, see Fig 30. 

15 

5.7 Cross-Hybridization Analysis: As described earlier ARROGANT estimates the amount of 
cross-hybridizing signal for each element on the array and tags the potential false positives. Fig 
31a-b shows the use of ARROGANT to design expression rnicroarrays. ARROGANT can be 
used to determine the regions that might cause cross hybridization and the researcher could 

20 design primers so as to exclude these regions which could have lead to false positives. 

ARROGANT can be used to tag the potential false positives due to cross hybridization. This 
helps the researcher in interpreting the expression data in a more thorough manner. 

5.8 Search using SQL: Reseachers can directly use various simple SQL queries to retrieve 

25 information from many different databases simultaneously. Also conversions like getting gene 
names for given accession numbers or accession numbers for given descriptions can be 
achieved using simple SQL queries. For example, using a single SQL query, LocusLink 
identifiers for all the 15,000 gene names in a collection can be achieved. 

30 Section 6: Results and Validation. 

ARROGANT was used for a variety of different applications, including cancer studies, 
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and clones were ordered based on the analysis done by ARROGANT. ARROGANT's ability 
to estimate potential cross-hybridization was tested on our human array containing 1 0,000 
elements and the cross-hybridizing candidates were categorized as either ALUs , simple 
sequence repeats or others. ARROGANT was used to recommend new genes for the Program 
5 for Genomic Application (PGA) and also used to annotate the existing genes on PGA. 
ARROGANT identified new genes for studies of Robert's Syndrome and annotated the 
associated gene list. The analysis mode of ARROGANT was used to annotate several different 
gene lists which included genes on commercial microarrays (Affymetrix, Inc.), and genes on 
our 10,000 human array and 12,000 mouse array. These and other applications using databases 
1 0 developed with ARROGANT are discussed below. 

6.1 ARROGANT used in cancer studies to recommend candidate genes, look for clones and 
determine the ones needed to be ordered: The utility of ARROGANT was demonstrated using 
all its three modes of operation to develop a new human microarray for cancer. Starting with 

15 the keyword search feature in the design mode, ARROGANT was used to compile a list of 
candidate genes for cancer. This list and two more lists obtained from separate sources were 
then combined to obtain a unique list using the merging gene collection mode. The genes 
already present on our human microarray (consisting of 10,000 genes) were then eliminated 
from the list. Now using the analysis mode, clones commercially available with Research 

20 Genetics were determined for the ones that needed to be ordered. ARROGANT was later used 
to verify the insert sizes of the clones obtained from Research Genetics. 

6.1.1 Keyword search to obtain possible candidate genes: The keyword search feature of 
ARROGANT was used to compile a collection of genes possibly associated with cancer using 

25 keywords like cancer, metastasis, carcinoma, etc. Collections of 359 probable candidate genes 
were selected. Now the Analysis mode was used to determine the UniGene clusters for this list 
of 359 genes to obtain a set of 237 non-redundant genes. 

6.1.2 Merging Gene Collection Mode: For a more complete coverage of candidate genes it 
30 was decided to merge the list obtained from ARROGANT with 2 more important collections 

from separate sources. These included 750 genes obtained from experts and another list of 
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13,969 genes from the CGAP website. The merging gene collection mode was used to compile 
one unique list of genes by avoiding sequence redundancies. The sequence redundancies were 
eliminated by having one gene per UniGene cluster. The merge list was compared to our 
human 10,000 microarray to determine the ones needed to be ordered. It was found that 9,315 
5 genes were not in common with our human 1 0,000 microarray. 

6.1.3 Look for commercially available clones with research genetics: Clones were found with 
research genetics for 7,593 out of 9,3 1 5 genes needed to be ordered. The clones included both 
sequence verified and non-verified. The results were found to be identical when compared to 

10 the output obtained from Research Genetics software program CMiner to look for 
commercially available clones. The 7,593 clones were ordered. 

6. 1 .4 Verification of clones obtained by determining their insert size: Once the clones were 
obtained from Research Genetics it was required to verify their insert sizes. The insert sizes 

15 can be determined only for the sequence verified clones as obtained from the Research 

Genetics database. There were 5,739 sequence verified clones out of the 7,593 clones received. 
ARROGANT could determine the insert sizes for 3,726 clones, as Research Genetics does not 
provide the insert sizes for all the sequence verified clones. 

20 6.2 Cross-hybridization analysis of the human array containing 10,000 elements. The results 
obtained from the experiments performed on yeast microarrays suggested that the threshold for 
a sequence to cross hybridize was an overlap consisting of 65 contiguous hydrogen bonds. 
ARROGANT was used to further study the effects of cross hybridization using human 
microarrays. Each gene on the 10,000 human microarray was BLASTed against the rest of the 

25 genes on the microarray. The BLAST results were used to find the potential 

cross-hybridization of each gene with the remaining 9,999 elements on the microarray. 
ARROGANT also determined the sequences having ALU elements in the 3' / 5' / coding 
region. Table 6.2 gives the distribution of the ALU elements. 



Type 


3* 


Coding 


5' 


ALU 


433 


0 


199 



Table 6.2: Distribution of ALU elements on our human 10,000 microarray. 



40 



UTSD:668 



There are a significant number of ALU elements in the 3' region. Most of the top candidate 
genes for cross-hybridization were found to contain ALU elements. The significant number of 
ALU elements in the 3' region suggests that there might be significant amount of 
cross-hybridization. ARROGANT calculates potential cross-hybridization of each element by 
5 sequences withn the chip as well as the entire UniGene database. 

6.3 ARROGANT used for identifying and annotating genes for polymorphism discovery to 
link to cardiac diseases for PGA: The Program for Genomic Application (PGA) is a 
nationwide attempt to use genomic and proteomic methods to study and investigate cellular 

10 responses to injury and inflammation. The program endeavors to identify the genes and 
proteins involved in these responses. ARROGANT was used to both recommend new 
candidate genes for PGA as well as annotate the current PGA list of 253 genes. The ability of 
ARROGANT to find potential candidates was tested by comparing the list obtained using 
keyword search with the current list of genes. The list of keywords compiled by researchers 

15 participating in PGA was as follows: 



hyperlipidemia 


arteriosclerosis 


low density lipoproteins 


cholesterol 


dietary responsiveness 


inflammation 


high density lipoproteins 


cytokine 


coronary calcification 


orphan receptor 


insulin resistance 


cardiac failure 


cardiac hypertrophy 


signal transduction 


coronary artery disease 


G-protein 


coronary atherosclerosis 





25 

ARROGANT found 3,789 genes associated with the above keywords. There were 13 genes 
found in common with the current PGA list of 253 genes. This demonstrated the keyword 
search capability of ARROGANT to look for potential candidates associated with keywords. 
The newly compiled list was annotated using the analysis mode and is available on the web at: 
30 http://ARROGANT.swmed.edu/myweb/hideandsort.asp?txt_array=40710. ARROGANT was 
also used to annotate the current PGA list of 253 genes. 
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The ability of ARROGANT in the analysis mode to accept a list of genes tab delimited 
with a number was used to assign priority levels to the genes: 2- Highest priority, 1- Moderate 
priority and 0- Low priority. The annotated table is available on the web at: 
http :// ARROGANT. swmed. edu/myweb/hideandsort. asp?txt_array=44082 . 

6.4 ARROGANT used in the study of Robert's Syndrome: Robert's Syndrome is a genetic 
disorder caused by chromosome damage during cell division, and characterized by loss of limb 
bones, cleft palate, heart defects and abnormalities of the abdominal organs. ARROGANT was 
used to find new potential candidate genes for Robert's syndrome using keywords: 
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Robert syndrome 


hypoplastic nasal and auricular cartilage 






atrial septal defect 




Robert's syndrome 


patent ductus arteriosus 




Pseudothalidomide syndrome 


polycystic kidneys 


5 


SC phocomelia syndrome 


fused kidneys 




bet ero chromatin 


horseshoe kidneys 




j^gterochroinatic repulsion 


micronucleation 






enlargement of the phallus 




Premature centromere separation 


absent nails 


10 


premature separation 


ICF syndrome 




Tetraphocomelia 


Centromeric instability immunodeficiency 




Limb reduction 


syndrome 




hypoplastic 


MECP2 




Long bone 


Methyl binding protein 


15 


Aneuploidy 


Hypometfiylation 




Craniofacial 






Oxycephalic 


Demethylation 




aplasia of the fibula 


demethyltransferase 






Methyl ation 


20 


a senceo 


methyl ase 




clett lip and palate 


mSIN3A 




oligodactyly 


Histone 




microcephaly 


Historic acetylation 




exophthalmus 


Histone acetylase 


25 


hypertelorism 


Histone deacetylase 




corneal clouding 


TAR syndrome 




hemangiomas 





ARROGANT found 6,326 genes, which were further annotated using the analysis mode. The 
30 results are available on the web at: 

http:// ARROGANT. swmed.edu/myweb/hideandsort.asp?txt_array=12345. A separate list of 
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16 gene names found to be important in the study of Robert Syndrome was obtained. The 
accession numbers for these 1 6 genes were determined using ARROGANT. It was found that 
there was one gene in common between the two lists. This again demonstrated the utility of 
ARROGANT to look for and identify candidate genes associated with keywords. The list of 16 
5 genes was also annotated using the analysis mode and the results are available on the web at 
http :// ARROGANT. swm ed.edu/myweb/hideandsort. asp?txt_array=3 96 1 3 . 

6.5 ARROGANT used to annotate genes on commercial DNA chips: ARROGANT was used 
in the analysis mode to annotate various microarrays available from Affymetrix (Santa Clara, 

1 0 C A) to help the researcher view the results obtained from the expression studies in a 
convenient manner. This provides the researcher a group of genes having particular 
characteristics together and helps in making important observations. The following 
commercial (Affymetrix) human and mouse microarrays were analyzed. 

1. Human HUG95 microarray: This microarray consists of 12,454 different elements. 
1 5 The annotated list is available on the web at 

http://ARROGANT.swmed.edu/myweb/hideandsort.asp?txt_array=l 1111. 

2. Rat RG-U34 microarray: This consists of 1,322 genes from Rat genome. The results 
are available on the web at 

http://ARROGANT.swmed.edu/myweb/hideandsort.asp?txt_array=57860. 

20 

6.6 ARROGANT used to annotate genes on chromosome 3p: ARROGANT was used to 
identify genes commonly mutated or whose expression is deregulated in human lung and 
breast cancers. Although several regions of loss occur on multiple chromosomes it was 
observed that allele loss in the chromosome 3p21.3 area was the earliest pre-malignant change 

25 so far detected in lung cancer development 

(http://www.utsouthwestern.edu/cancer/Research/3p21_intro.htm). ARROGANT was used to 
annotate the 32 genes on chromosome 3p thought to be important in causing lung cancer. The 
results are available at: 

http : / /ARROGANT . s wmed . edu/myweb/hideands ort . asp ?tx t_array=403 5 7 . 

30 

6.7 ARROGANT used to analyze human microarrays: Our laboratory has developed a human 



44 



UTSD:668 



cDNA microarray, which consists of 10,000 clones from Research Genetics. Many 
laboratories in UTSW (University of Texas Southwestern Medical Center at Dallas) are using 
this microarray for various research studies like cancer, aging, etc. ARROGANT provides 
annotation for all the genes as one table. The researchers can overlay their expression level 
5 data on this table, which would help them make important observations. For example, the 

researcher could look at the pathways for all the highly expressed genes and also know their 
position in the genome. Further the researcher could also sort the data using ARROGANT to 
bring the interesting genes on top of the table. ARROGANT annotation of the human 10,000 
array is available on the web at 
10 http://ARROGANT.swmed.edu/myweb/hideandsort.asp?txt_array=601 10. ARROGANT also 
annotated our earlier human array consisting of 4,200 elements and the results are available at 
http :// ARROGANT . swmed. edu/myweb/hideandsort. asp?txt_array=607 1 8 . 

6.8 Other Applications developed writing simple SQL queries: Due to the implementation of 
15 various databases locally in a relational database format, certain features from across databases 
can be easily retrieved by writing simple SQL queries. ARROGANT has been used for many 
such applications including: 

1 . Genomic entries and LocusLink: identifiers for 793 gene names were required. 
LocusLink identifiers were first retrieved for each of the gene names and genomic entries for 

20 each of the LocusLink identifiers were retrieved. This was possible by writing single SQL 
query since all the databases were implemented locally and stored in the relational database 
format. 

2. Another application required finding pathways for 521 genes having GenBank Numeric 
Identifiers (NfD). Pathway information for each of the 521 genes was retrieved by writing an 

25 SQL script using ARROGANT database. 

From the above discussion it is seen that ARROGANT can be used in a variety of different 
applications ranging from annotating any list of genes, recommending new potential 
candidates associated with keywords to merging different lists of genes to obtain one 
collection without redundancies. 

30 

Section 7: Enhancements and add-ons. 
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ARROGANT has been applied to a variety of different applications as demonstrated in 
the earlier section and has proved to be a worthy tool. It compiles, annotates and merges large 
gene collections as well as helps in the design of expression / resequencing microarrays. The 
following optional features and add-ons enhance its ability to work with large gene collections 
5 and improve its ability as a design and analysis tool for microarrays. 

1. Include more databases: ARROGANT uses multiple databases in a modular format 
and makes it convenient to add more databases. Parsers may be used to add more databases 
like Swissprot and Pfam, to increase the annotating capability of ARROGANT for large gene 
collections. Additional clone databases, such as proprietary databases (e.g. LifeSeq, Incyte 

10 Genomics) enhance the ability of ARROGANT to look for commercial clones, and the 

addition of literature databases like Medline (NCBI) increase the value of ARROGANT to 
look for potential candidates by keyword search. 

2. Ability to find introns/exon boundaries: ARROGANT can implement programs like 
GENSCAN (Burge et al, J. Mol. Biol. 1997) to detect introns / exon boundaries for a given 

15 gene sequence, enabling the tool to detect possible splice variants and design primers to 
amplify each exon. 

3. Clustering based on pathways: Basic ARROGANT has the capability of sorting the 
final annotated table (analysis mode) based on pathways. This program is enhanced to use 
pathway information to cluster together up-regulating / down-regulating elements on the array 

20 for each individual pathway. 

4. Include keyword search within analysis mode: The user could use the analysis mode 
to annotate a large gene collection and then use the keyword search within the analysis mode 
to display only those genes related to a certain keyword. For example, the researcher would 
want to know the genes associated (by keywords) with cancer from the ones over-expressed in 

25 a given microarray experiment. This could be achieved by letting the researcher type in the 

keywords to be searched and writing a function to look for those keywords within the current 
annotation table of the analysis mode. 

5. Integrate ARROGANT with other processes: Basic ARROGANT uses other 
laboratory tools like cross-hybridization (to estimate amount of false signal obtained due to 

30 genes having significant sequence similarity) and Rep-X (to find potentially polymorphic 

repeats). The output of more tools like SNIDE which may be used to predict candidate SNPs 
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(Single Nucleotide Polymorphism), may be added as a column to the annotation table of the 
analysis mode, i.e. adding the output of SNIDE as a table in the ARROGANT database. 

6. Use of English query language: SQL Server 7.0 can be trained to enable the users 
(having direct access to the database) to type queries in simple English language instead of 

5 SQL (Structured Query Language). For example SQL server 7.0 software can be trained so 

that the user could type in queries like ' Give me genes related to cancer but not aging' instead 
of writing a query programmatically in SQL 

7. Include expression level data for more than one experiment in the analysis mode: 
ARROGANT can also be made to provide basic statistical results e.g. mean, max, etc. along 

10 with other annotation. This is to say that expression data from more than one experiment may 
be included in the same annotation table and columns giving the mean, max and median can 
be added. In this case ARROGANT directly inputs gene collection, tab delimited by one or 
more expression level data, and the program automatically determines number of experiments 
by counting the tabs between each expression data provided. 
15 8. Enhance keyword search: Basic ARROGANT keyword search can be used with two 

r. Boolean operators AND/OR'. The keyword search capability can be further extended to 

include more boolean operators like NOT, BUT etc. e.g. Aging NOT telomerase' which would 
^ translate to look for genes associated with aging but not with telomerase. 

9. Basic ARROGANT accepts accession numbers as input identifiers; enhancements 
20; provide for accepting gene sequences using various identifiers e.g. GenBank NID, UniGcne 
cluster ID, gene names and LocusLink ID. 

All publications and patent applications cited in this specification and all references 
cited therein are herein incorporated by reference as if each individual publication or patent 

25 application or reference were specifically and individually indicated to be incorporated by 
reference. Although the foregoing invention has been described in some detail by way of 
illustration and example for purposes of clarity of understanding, it will be readily apparent to 
those of ordinary skill in the art in light of the teachings of this invention that certain changes 
and modifications may be made thereto without departing from the spirit or scope of the 

30 appended claims. 
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